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For a model deca-alanine peptide the cavity (ideal hydrophobic) contribution to hydration favors 
the helix state in the coil-to-helix transition and the paired helix bundle in the assembly of two 
helices. The energetic contributions of attractive protein-solvent interactions are separated into a 
short-range part arising from interactions with solvent in the first hydration shell and the remain¬ 
ing long-range part. In the helix-coil transition, short-range attractive protein-solvent interactions 
outweigh hydrophobic hydration and favor the unfolded coil states. Analysis of enthalpic effects 
shows that it is the favorable hydration of the peptide backbone that favors the unfolded state. 

Protein intramolecular interactions favor the helix state and are decisive in folding. In the pairing 
of two helices, the cavity contribution outweighs short-range attractive protein-water interactions. 

However, long-range, protein-solvent attractive interactions can either enhance or reverse this trend 
depending on the mutual orientation of the helices. In helix-helix assembly, change in enthalpy 
arising from change in attractive protein-solvent interactions favors disassembly. In helix pairing as 
well, favorable protein intramolecular interactions are found to be as important as hydration effects. 
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Helices have been termed the “hydrogen atoms of 
biomolecular complexity” [l[ . In proteins the a-helix is a 
common structural motif and understanding the forma¬ 
tion of a-helices occupies a pre-eminent place in efforts 
to understand protein folding. Using computer simu¬ 
lations and a new approach to free energy calculations 
0, i, here we revisit this classic problem and study 
two transitions in a model deca-alanine peptide. Mirror¬ 
ing the primary-to-secondary and secondary-to-tertiary 
transitions in protein folding, we study, respectively, the 
coil-to-helix transition and the pairing of helices to form 
a helix dimer. Our principal focus is to explicate the 
hydrophobic and hydrophilic hydration contributions in 
these transitions at a level that has hitherto been possible 
only for small molecular solutes. 

There are several reasons to re-examine the fundamen¬ 
tal premises about hydration thermodynamics of pro¬ 
teins and the forces driving protein folding. First, re¬ 
cent experiments and simulations challenge the conven¬ 
tional view that hydrophobic interactions drive protein 
folding. These studies show that the all-backbone polyg¬ 
lycine, and analogous archetypes of intrinsically disor¬ 
dered peptides, can undergo a collapse transition in wa¬ 
ter UQ- Second, recent computer simulations 0-[l3 
reveal important physical and conceptual limitations in 
the group-additive approach that has been the back¬ 
bone of approaches to understand protein hydration ther- 


* Corresponding author: Dilip.Asthagiri@rice.edu 


modynamics in experiments (cf. Ref. [T^ and [ii). In¬ 
deed the prevailing views of protein hydration thermo¬ 
dynamics [iM3 trace back to this group-additive rea¬ 
soning. Lastly, while simulations can in principle pro¬ 
vide a detailed molecular thermodynamics understand¬ 
ing of hydration, with some exceptions [3 m , such 
studies are scarce for realistic proteins and polypeptides. 
The availability of an approach that alleviates this situ¬ 
ation presents an opportunity to revisit and potentially 
re-appreciate a classic problem in protein folding. 

Earlier studies based on continuum solvent or lattice 
models have come to differing conclusions about solvent 
effects in the coil-to-helix transition. Some have sug - 
gested that hydrophobicity drives the transition [3 lUl , 
while others have emphasized the role of favorable elec¬ 
trostatics [12. Experiments suggest that helix extension 
is enthalpically driven [23l | , as has also been found in com¬ 
puter simulations interpreted within the Zimm- Bragg o r 
Lifson-Roig formalisms (for example, see Refs. I24l - |2CT . 
Interestingly, both experiments |27| and computer simu¬ 
lations |2J-l26j show a negative heat capacity upon un¬ 
folding, the opposite of what is observed in unfolding of 
globular proteins [12. The negative heat capacity has 
been interpreted as arising due to the favorable hydra¬ 
tion of the backbone upon unfolding [12 , but a molecular 
scale description of this signature remained to be sought. 

The approach we have developed is based on a qua¬ 
sichemical organization [lol - llTj l of the potential distri¬ 
bution theorem [12 . Using this approach, we are able 
to interrogate the hydration thermodynamics of proteins 
ii at a level of resolution that is comparable to those 
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for small molecular solutes. Our study on cytochrome 
C helped reveal limitations of continuum models of hy¬ 
dration Q. Subsequent studies have helped illuminate 
conceptual and physical limitations in the group-additive 
description of the hydration thermodynamics of the pep¬ 
tide backbonefllj and of a hydrophobic side-chain in the 
context of model peptides [T^l ■ On the basis of these new 
developments, and for the reasons noted above, here we 
study the hydration thermodynamics in the coil-to-helix 
and helix-helix assembly in a deca-alanine peptide. 

We find that hydrophobic interaction is not consequen¬ 
tial in the formation of the helix, but it does play an 
important role in helix-dimerization. But hydrophilic 
hydration is found to play a nontrivial role in the coil- 
to-helix transition and helix-dimerization. Indeed hy¬ 
drophilic hydration can drive the unfolding of the helix 
and disassembly of the helix dimer. Throughout back¬ 
bone hydration is found to be the most important com¬ 
ponent of the overall hydrophilic hydration. For the sys¬ 
tems studied here, hydrophilic hydration and protein in¬ 
tramolecular interactions are as important as, if not more 
important than, hydrophobic effects. 


I. THEORY 

The excess chemical potential, is that part of the 
Gibbs free energy that arises from intermolecular inter¬ 
actions and of first interest in understanding the solubil¬ 
ity of a solute in a solvent. Here /r®’' is defined relative 
to the ideal gas at the same density and temperature. 
To calculate /i®’', we regularize the statistical problem 
of calculating it from the potential distribution theorem 
= ln(e^®) [13, 113 , where the averaging {...) is over 
the solute-solvent binding energy (e) distribution P{e). 
As usual P = l/feeT, with T the temperature and the 
Boltzmann constant. 

We introduce an auxiliary field ^(A; r) that moves the 
solvent away from the solute to a distance A. The dis¬ 
tance between the center of the field and the solvent 
molecule is r. For r > A, (/) = 0. Since the solvent inter¬ 
face is pushed away from the solute, the solute-solvent 
interaction is tempered and the conditional distribution 
P(£|(/)) is better behaved than P{£) iiini- In practice, 
we adjust the range A such that F’(£|(/)) is a Gaussian. 
With the introduction of the field, we havefi, H, [Tl| 

/3/r®’'= + pfi^^[P{£\(l))]- lnpo[(j)] ^ (1) 

local chemistry long-range packing 

where —fcB’rina;o[(/)(A)] is the work done to apply the 
field in the presence of the solute, —fcBrinpo[</'(A)] is the 
corresponding quantity in the absence of the solute, and 
/3^®’'[P(£|^)] is the contribution to the interaction free 
energy in the presence of the field. Fig. [T] provides a 
schematic description of Eq. [T] 

We apply the field about each heavy atom to carve 
out the molecular shape in the liquid (Fig. [1}. For con¬ 
venience we use the same value of A for all the heavy 
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FIG. 1. Quasichemical organization of the excess chemical 
potential. The A = 3 A envelope defines the solvent excluded 
(SE) volume and the A = 5 A defines the envelope extending 
to the first hydration shell (FS). For chemistry coupled (un¬ 
coupled), the solute-solvent interaction is present (absent). In 
Eq. [T] we follow FS-packing to the hydrated solute; in Eq. [5] 
we follow SE-packing. Figure adapted from Ref. [TJ with per¬ 
mission from Elsevier. 


atoms. For A > 5 A, P{£\(j)) is well-described by a Gaus¬ 
sian. We thus define Ag = 5 A. 

The solute also excludes a volume to the solvent and by 
definition there are no short-range attractive interactions 
between the solute and the solvent in this domain. We 
find that Inio ~ 0 for A < 3.0 A, irrespective of the con¬ 
formation of the peptide. This suggests that the space 
enclosed by Ase = 3.0 A is excluded to the solvent and 
thereby provides a natural definition of molecular extent 
of the cavity to be used in discussions of the cavity (ideal 
hydrophobic) contribution. Interestingly, the range be¬ 
tween 3 A to 5 A corresponds to the first hydration shell 
for a methyl carbon [s^l and is a conservative description 
of the hrst hydration shell of groups containing nitrogen 
and oxygen heavy atoms. For simplicity 

Using the Ase and Ag, we rearrange Eq. [T]as 


/3/r®’' = In 


a;o[0(AG)] 

Po[</'(Ag)]/po[<('(Ase)] 


revised chemistry 


-lnpo[</’(AsE)] ■ 

'-V-' 

SE packing 


/3/.®=^[P(£|</>(Ag))] 

'-V-' 

long-range 

( 2 ) 


Note that SE packing and the revised chemistry plus 
long-range contribution is uniquely defined for the given 
forcefield. Physically, the revised chemistry contribution 
measures the free energy contribution from solute inter- 
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action with solvent in the first shell relative to a non¬ 
interacting (uncoupled) solute (Fig[T|). 

The revised chemistry and long-range contributions de¬ 
scribe the role of short-range and long-range attractive 
protein-solvent interactions in the thermodynamics of hy¬ 
dration. These two components of the hydrophilic contri¬ 
bution occur at different energy (and length) scales. But 
within commonly used continuum models of electrostatic 
interactions this distinction is necessarily lost and both 
the short-range and long-range interactions are treated 
as part of long-ranged interactions (e.g. see da). From 
the perspective of such continuum models, our definition 
of long-range interaction is more conservative. 

The SE packing contribution measures the hydropho¬ 
bic hydration of an ideal hydrophobe [nisi.' In theo¬ 
retical discussions of hydrophobic effects, a cavity with 
a hard-wall interaction in water is often considered. The 
packing contribution in our calculation uses a soft-cavity. 
This soft-cavity packing estimate is always a lower-bound 
to the hard-cavity estimate and can be easily corrected 
to give the latter y . We do not pursue those corrections 
here and instead use the soft-cavity packing result as a 
measure of primitive hydrophobic effects. We refer the 
reader to published papers iiini for more extensive 
details about the approach. 


A. Helix-Helix PMF 

The above development carries over to the calculation 
of kF(r), the potential of mean force (PMF) to bring 
two helices a distance r apart, where r is the separation 
between the helix axis with the axis parallel to each other. 
The PMF 

W{r) = Wsoiv{r) + ^U[r), (3) 

where Wsoiv is the solvent (or indirect) contribution 
[ 3 ^ and A {7 is the contribution from direct protein- 
protein interactions. Wsolver) = Md?mer(^) “ 2^®^, where 
A^dtmer(^) hydration free energy of the pair of he¬ 

lices (for a given separation and orientation) and /i®^ is 
the hydration free energy of a monomer helix. 


B. Entropic and enthalpic contributions 


For calculating Ereorg, we adapted the hydration-shell- 
wise calculation described earlier [H,[sj (cf. Sec. S.III). 


H. METHODS 


The regularization follows previous work m- The 
deca-alanine peptide was modeled with an acetylated 
(ACE) N-terminus and n-methyl-amide (NME) capped 
C-terminus. The extended /3-conformation (cj), ip = 
— 154 ± 12,149 ± 9) was aligned such that the end-to-end 
vector lay along the diagonal of the simulation cell. We 
label this coil state as Cq. The helix was aligned with 
the helix-axis along the x-axis of the cell. The initial 
structures were energy minimized with weak restraints 
on the heavy atoms to relieve any strain in the structure. 
The peptides were solvated in 3500 TIP3P 36,[37| water 
molecules. Version c31 of the CHARMM 381 forcefield 
with correction(cmap) terms for dihedral angles [s^, was 
used for the peptides. 

We sampled unfolded states using the adaptive-bias 
force (ABF ( 4 ^, 1 ^) approach which also additionally pro¬ 
vided the free energy of unfolding the polypeptide in vac¬ 
uum. From the ABF trajectory, we sampled nine struc¬ 
tures with end-to-end distances between terminal carbon 
atoms ranging between 28 A and 36 A in increments of 
1 A. We label the coil states from this unfolding simu¬ 
lation {Cl,..., Cg}. The (j), Ip for these unfolded states 
predominantly populate /3 and PPII regions of the Ra- 
machandran plot. (Note that, in the strict sense of coil- 
to-helix transition, the extended /3 is one extreme mem¬ 
ber in the ensemble of coil states.) In all the hydration 
free energy calculations, the structures were held rigid. 

The same set-up was used to investigate helix pairing. 
Additionally, we consider two relative orientations of the 
helix dipoles, parallel and antiparallel. (Note that the he¬ 
lix dipoles will be antiparallel in the simplest helix-turn- 
helix motif.) These arrangements help illuminate the role 
of long-range protein-solvent interactions in helix-helix 
complexation. We note that in nature helices almost 
never align perfectly parallel or antiparallel [i^ , but the 
hydration effects that are of principal interest here are 
expected to be insensitive to minor distortions in the rel¬ 
ative orientation. 


Ignoring the pressure-volume correction and contribu¬ 
tion due to a finite isothermal compressibility of water, 
the enthalpy of hydration, /i®’', is given by 

/l®’' = Esw + Ereorg (4) 

where Esw is the average peptide-solvent interaction en¬ 
ergy and Ereorg is the (average) water reorganization en¬ 
ergy. Additionally ignoring the contribution due to a 
finite thermal expansivity of water, the entropy of hy¬ 
dration is given by 

( 5 ) 


HI. RESULTS AND DISCUSSION 
A. Coil-to-helix transition 

Figure [2] collects the results of the hydration free ener¬ 
gies of the helix and (Cq, ..., Cg} coil states; it is clear 
that the coil states are better hydrated than the helix. 
In Table |T] we present the free energy components of the 
helix state, the least favorably hydrated coil state (Cq), 
and the most favorably hydrated coil state (C 7 ). 

Before we discuss our results, we first compare our cal¬ 
culations with existing hydration free energy data. As 


Ts®’' = . 
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TABLE 1. Components of the hydration free energy for the helix and the least favorably (Co) and most favorably (C 7 ) hydrated 
coil states. For the helix and Co states, results with partial charges turned off (indicated by Q = 0) are also noted. Rg is 
the radius of gyration (relative to the center of mass) and Rc is the end-to-end distance between terminal carbon atoms in A. 
SASA is the solvent accessible surface area in A^. All energy values are in kcal/mol. Standard error of the mean is given in 
parenthesis at the 2 ct level. 


Conformation 

Rg 

Rc 

SASA 

SE 

Revised Chemistry 

Long-range 

p- 

Ereorg 

Esw 


Helix 

5.3 

16.3 

876.0 

44.4 (0.4) 

-51.6 (0.2) 

-31.6 (0.2) 

-38.8 (0.5) 

75 ( 6 ) 

-150.9 (0.6) 

-37.1 

Helix (Q = 0) 





-22.4 (0.3) 

-18.1 (0.03) 

3.9 (0.5) 

42 (4) 

-59.7 (0.6) 

- 21.6 

Coil (Co) 

11.1 

36.8 

1260.0 

58.4 (0.3) 

-77.1 (0.3) 

-27.6 (0.1) 

-46.3 (0.4) 

92 ( 6 ) 

-186.8 ( 0 . 6 ) 

-48.9 

Coil (Q = 0) 





-30.9 (0.4) 

-24.7 (0.03) 

2.8 (0.5) 

50 ( 8 ) 

-80.8 ( 0 . 6 ) 

-33.6 

Coil (C 7 ) 

10.6 

34.0 

1249.0 

58.3 (0.5) 

-83.4 (0.3) 

-28.3 (0.4) 

-53.4 (0.8) 

98 (4) 

-201.4 (0.6) 

-57.4 



FIG. 2. Hydration free energies of the helix and {Co,..., Cg} 
coil states. The horizontal axis has no meaning and is used 
solely to differentiate multiple coil states with similar val¬ 
ues. The radius of the symbol is equal to twice the standard 
error of the mean {2a). The Co state is the smallest circle in 
the collection of unfolded states; the standard error is about 
half compared to the other estimates because we had 4 times 
more data for Co (Sec. S.I). 


already noted above, till recently very few all-atom cal¬ 
culations of the hydration free energy of polypeptides 
with around 10 residues appears to have been reported. 
Helms and coworkers (isj have studied blocked-(Ala)„ he¬ 
lix hydration with n up to 9. Extrapolating their results 
suggests a value of about —40 kcal/mol for the deca- 
alanine helix, consistent with the quasichemical results 
(Table HI). Kokubo et al.[l^, computed the van der Waals 
and electrostatic contribution to the hydration of a deca- 
alanine helix and for several coil states. Our hydration 
free energies based on their structures are in fair agree¬ 
ment (data not shown), and the agreement becomes ex¬ 
cellent with more extensive sampling in the calculation 
of the van der Waals contribution [^. Besides these, 
the regularization approach has been tested in studies 
on water [ 43 . l45l|. ions [4^, and simple peptides [ll|, and 
the protein cytochrome ci. Further, the quasichemical 
framework has also been thoroughly documented (29l - l3ll| . 


Turning to the results of this study. Table |T] indicates 
that at least for the coil states considered here, the pack¬ 
ing is somewhat insensitive to the peptide structure. At 
the scale of the cavities, we expect t he p acking contribu¬ 
tions to scale with surface area [43-1^ and the data in 
the table conforms to this expectation. The free energy 
per unit SASA is similar to values rep orted for hard- 
spheres using scaled-particle theory [43, and explicit 
all-atom simulations [Sll - I^ , but is about a factor of ten 
larger than the surface energy parameter often used in 
continuum surface-area based model of hydrophobic hy¬ 
dration [ 5 ^ . Since the solvent-accessible surface areas of 
Co and C 7 are not very different, despite overall differ¬ 
ences in the structure, we expect the packing contribu¬ 
tion to be similar for these states, as found in simulations. 
The packing contribution favors the helix state by about 
— 14 kcal/mol: as expected, hydrophobic hydration fa¬ 
vors the compact state of the protein. 

The chemistry contribution, however, favors the coil 
states by between —31.8 kcal/mol for the C 7 state and 
—25.5 kcal/mol for the Co state. Thus, the local protein- 
solvent interaction outweighs the packing contribution 
by between 11 kcal/mol to 18 kcal/mol in favoring the 
coil state. Comparing the chemistry contribution for the 
peptide and its Q = 0 analog (Table |T| shows that the 
favorable local protein solvent interactions arise primar¬ 
ily from favorable peptide backbone-water interactions, 
the role of the methyl groups in electrostatic interactions 
being comparatively negligible. This observation is di¬ 
rectly confirmed from the analysis of the solute-solvent 
interaction contribution to the enthalpy (Table IJ). 

TableUalso shows that the entropy of hydration is neg¬ 
ative, but here it arises due to attractive solute-solvent in¬ 
teractions. Anticipating our forthcoming article on tem¬ 
perature effects, we note that the entropy calculated us¬ 
ing Eq. [ 5 ] is in agreement within statistical uncertainties 
with calculated from the temperature derivative of 
as is expected for thermodynamic consistency. 

Results for the Q = 0 analog suggests that attrac¬ 
tive solute-water dispersion interactions alone can inhibit 
chain collapse brought about by packing effects (Table|I|. 
The enthalpy of hydration is large negative, despite the 
positive contribution from solvent reorganization (Ta¬ 
ble [J), emphasizing the role of attractive dispersion inter¬ 
actions between the Q = 0 peptide and solvent. These 
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observations emphasize that care is needed in assuming 
the relevance of the poor solubility of nonpolar solutes 
in rationalizing the collapse of a polypeptide. Our obser¬ 
vation that hydration does not explain the collapse of a 
nonpolar chain is consistent with the observation of sim¬ 
ilar behavior in alkanes (cf. Ref. and the reanalysis of 
data in Ref. [H presented therein). Interestingly, attrac¬ 
tive solute-water interactions also oppose the pairing of 
the prototypical hydrophobe methane [s^ : emphasizing 
the importance of solute-water attractions, this effects is 
more pronounced for larger alkanes [s^ . 

Protein-solvent long-range interactions (cf. Fig.[T|) con¬ 
tribute a substantial fraction — between 50% to 82% — 
of the net hydration free energy of the peptide. About 
90% of the (favorable) long-range contribution for the 
Co and C 7 states arises from van der Waals interactions, 
while for the helix it is about 60%. As can be expected, 
the electrostatic contribution for the helix is higher be¬ 
cause of the high dipole moment of the helix. (The long- 
range contributions are non-negligible because the num¬ 
ber of solvent molecules in the outer-shell domain is not 
small, although individual water-solute interaction ener¬ 
gies are typically small.) The long-range electrostatic and 
van der Waals contributions balance in the coil-to-helix 
transition resulting in a net free energy change of about 
—4.0 kcal/mol in favor of the helix (Table U). 


The above analysis shows that the hydrophilic contri¬ 
butions to hydration outweigh the hydrophobic driving 
force and favor the unfolded state of the peptide. Exper¬ 
iments show that a coil-to-helix transition in a predomi¬ 
nantly alanine-based peptide can occur for a polypeptide 
chain comprising as few as 13 residues [l^. We expect 
the role of hydration in disfavoring the coil-to-helix tran¬ 
sition to hold for this slightly larger chain, provided the 
coil states are such that the backbone remains accessi¬ 
ble to solvent. This then suggests that the experimen¬ 
tally observed coil-to-helix transition must be driven by 
changes in protein intra-molecular interactions, an in¬ 
ference that is in consonance with the suggestion that 
enthalpy changes driving helix formation 23|, albeit in 
longer chains. Results in Table HIl support this expecta¬ 
tion. 


Table in shows that the favorable hydration of the 
backbone is lost in coil-to-helix transition {AEbb > 0 ) 
and this is larger than the change in enthalpy of solvent 
reorganization {AEreorg < 0 ) which favors the more com¬ 
pact helical state. It is clear that a favorable change in 
the internal energy {AEint < 0 ) is necessary to obtain a 
favorable change in the net enthalpy Ah™ < 0 . 

On a per-residue basis the net change in enthalpy 
in the coil-to-helix transition is estimated to be be¬ 
tween —2.4 ± 0.8 kcal/mol/residue (Cq) and —2.8 ± 
0.6 kcal/mol/residue (C 7 ), with statistical uncertainty 
reported at the 2a level. A direct comparison of our es¬ 
timated enthalpy change per residue with experiments is 
hampered by (a) the lack of a rigorous conformational 
averaging in our calculations, (b) the short length of 
our polypeptide, and (c) the fact that in experiments 


TABLE II. Components of the enthalpy change in the coil-to- 
helix transition for the coil states noted below. The change in 
excess enthalpy (Ah™) reveals the role of hydration and the 
change in internal energy (AEint) the role of intra-molecular 
interactions. Ah™, is further divided into AEreorg, the 
change in the water reorganization contribution, and AEsw, 
the peptide water interaction contribution. The latter is sub¬ 
divided into contributions from the backbone-water, AEhb, 
and sidechain-water, AEsc, interactions. Standard error of 
the mean is given in parenthesis at the Icr level. All values 
are in kcal/mol. 


Coil 

A F, 

i—^j^reorg 

AEhb 

^Esc 

AEint 

^^total 

Co 

C 7 

-17.0(8.0) 

-23.0(6.0) 

34.0(0.4) 

51.1(0.2) 

1.9(0.2) 

-0.6(0.2) 

-43.3 

-55.0 

-24.0(8.0) 

-28.0(6.0) 


there are residues besides alanine to aid in solubilizing 
the peptide. Experiments on predominantly alanine pep¬ 
tide (with about 50 residues) suggests a value of 
— I.O kcal/mol/residue. Theory suggests that coil states 
that are less well-hydrated, similar to the Cq state, will 
dominate the net hydration thermodynamics (cf. Ref. 
and [13) for a corresponding result for ions). These are 
also the states that appear to have an enthalpy change 
closer to the experimental result (within statistical un¬ 
certainties of the calculation). Eorceheld bias can be an 
issue but using the recently re-optimized variant 
(C36 [l^]) of the forcefield changes our results by only 
an additional 5%. 


B. Helix-helix complexation 

We next consider the free energy of helix association 
or the potential of mean force between two helices. Fig. [3] 
shows that hydration opposes the complexation of helices 
in the antiparallel orientation. Given that the repulsion 
starts at a considerable interhelical distance, our choice 
of helical registration is probably a minor effect. Our 
results suggest that hydration will oppose formation of 
the helix-turn-helix motif. Interestingly, the direct in¬ 
tramolecular contributions (AU) can outweigh the hy¬ 
dration effects to drive complexation. In the antiparallel 
arrangement favorable AU drives complexation, whereas 
for the parallel arrangement, unfavorable AU tempers 
the favorable hydration effects. Reminiscent of protein 
folding free energies, the net free energy of complexation 
is roughly —2 kcal/mol, a small magnitude relative to 
the large competing hydration and inter-molecular inter¬ 
action effects. 

Figure [4] shows that primitive hydrophobic effects do 
drive helix-helix complexation and, in contrast to the 
coil-to-helix transition, they do outweigh the local chem¬ 
istry contributions. Thus packing (hydrophobic) effects 
do become important at larger length-scales. However, 
the long-range hydration interactions play an entirely 
nontrivial role in the complexation process, despite the 
peptides being net neutral; these effects can easily out- 
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r (A) 


FIG. 3. Components of the potential of mean force in bringing 
two helices together. The helices are shown in green and the 
red arrows indicate their mutual orientation. Wsoiv is the 
solvent contribution (open symbols), and Wsoiv+ AU (Eq.|3]) 
is the net PMF (filled symbols). For r < 8 A, there is steric 
overlap between the helices and AU rises rather sharply. Data 
including these values of AU are thus not shown. 


weigh the hydrophobic contribution. The antiparallel ar¬ 
rangement of the helices is strongly disfavored by loss of 
favorable solute-solvent interactions in hydration, while 
the inhibition is more modest for the parallel arrange¬ 
ment of helices. 

Fig. [5] shows that the orientation dependence of the 
long-range contribution arises solely due to electrostatic 
interactions. While van der Waals interactions between 
the solute and solvent (outside in the first hydration 
shell) do not discriminate between the two orientations, 
its magnitude is non-negligible on the scale of the helix 
complexation free energy. The parallel arrangement of 
the helices is favored because of the synergistic effect of 
the helix macro-dipole; the loss of hydration of the helix 
dipole also explains the unfavorable contribution for the 
anti-parallel arrangement. 

Analysis of enthalpic and entropic effects in pairing 
shows that for both parallel and antiparallel configura¬ 
tions entropic effects favor complexation (Table Hill) , but 
enthalpic effects do not. However, the characteristics of 
the change in reorganization and interaction components 
emphasizes the need for caution in interpreting the en¬ 
tropic driving force in terms of changes in water struc¬ 
ture. For the antiparallel configuration, the change in the 
water reorganization energy favors helix pairing, but its 
effect is negligible for pairing of parallel helices. For both 
orientations, the loss of hydrophilic backbone-solvent and 
sidechain-solvent interactions inhibits helix association. 
Emphasizing the importance of electrostatic interactions 
between the backbone and the solvent, the backbone- 
solvent contributions are sensitive to the orientation of 
the helices, but the sidechain-solvent contributions are 


essentially of similar magnitude. 

IV. CONCLUSIONS 

We have considered the hydration contributions to 
both secondary and tertiary protein structure formation 
by considering the change of a solvent-exposed coil to 
helix and the association of two such helices in a model 
decaalanine peptide. The latter idealized model does ig¬ 
nore the role the loops connecting the helices play in the 
pairing, but it is a reasonable starting point to under¬ 
stand solvent effects in tertiary structure formation. 

The association of the cavities, the prototypical hy¬ 
drophobic interaction, tends to favor the compact state of 
the polypetide and favor both the coil-to-helix transition 
and helix-helix complexation. But in the coil-to-helix 
transition, hydrophilic effects (protein-water attractive 
interactions) easily overwhelm the hydrophobic contri¬ 
bution and favor unfolding of the peptide. Even for a 
discharged peptide, essentially a nonpolar chain, attrac¬ 
tive solute-solvent dispersion interactions suffice to favor 
the unfolded state. 

In the pairing of helices, hydrophobic interactions out 
weigh the short-range peptide-water interactions in favor¬ 
ing helix complexation. This occurs at a larger length- 
scale than the coil-to-helix transition of a single peptide. 
However, long-range protein-solvent attractive interac¬ 
tions, especially for the antiparallel arrangement of he¬ 
lices, outweighs the net effect of the packing plus short- 
range attraction contributions to favor the disassembly 
of the helices. 

In both coil-to-helix transition and helix-pairing, the 
predominant hydrophilic effects (in our model system) 
arise from the interaction of the backbone with water. 
This observed importance of the backbone appears to 
be consistent with recent studies that encourage a re¬ 
appreciation of the role of the backbone in protein folding 
(for example, see [ISO and lbSl) . 

We find in both the coil-to-helix transition and in 
the pairing of the helices in the antiparallel orientation, 
changes in the intra-molecular energy of the protein are 
essential in shifting the balance to the folded state. The 
limitations of the models and forcefield notwithstanding, 
our study suggests that in protein folding hydrophilic ef¬ 
fects and protein intra-molecular interactions are as im¬ 
portant as, if not more important than, hydrophobic ef¬ 
fects. 
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FIG. 4. The SE packing, revised chemistry, and long-range (Eq. [2)) contributions to the free energy of helix-helix complexation. 
The data is presented relative to two helices infinitely apart. At contact (r « 9.5 A) the packing effects outweigh the local 
chemistry effects, but this trend can easily reversed by long-range interactions. 


TABLE III. Enthalpic (h®’') and entropic contributions (TAs®’') to the change in the hydration contribution to the potential of 
mean force, Waoiv, as the helices are brought from 14.5 A to 9.5 A. The change in hydration enthalpy, A/i®’‘, is further divided 
into AEreorg, the change in the water reorganization contribution, and AE^w, the peptide water interaction contribution. The 
latter is subdivided into contributions from the backbone-water interactions, AEhb, and sidechain-water interactions, AEsc- 
AX = A(9.5 A) — A(14.5 A), for all X. Standard error of the mean is given in parenthesis at the Icr level. The value of 
rs®’‘(r = 14.5 A) is explicitly noted to contrast with the sign of AEreorg- All units are kcal/mol. 


Gonformation 

AEbb 

AEsc 

AEreorg 

A/i®" 

AWsotv 

TAs®" 

Ts®=‘(r = 14.5 A) 

Antiparallel 

24.7 (0.5) 

6.2 (0.1) 

-7.5 (5.0) 

23.4 (5.0) 

13.8 (1.2) 

9.6 

-70.2 

Parallel 

4.3 (0.4) 

4.1 (0.1) 

0.8 (5.0) 

9.2 (5.0) 

-0.1 (1.2) 

9.2 

-69.3 
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FIG. 5. The van der Waals (solid line) and electrostatic 
(open symbols, dashed lines) contributions to the free en¬ 
ergy. The long-range contribution to free energy is uniquely 
decomposable [HI . 1^ into electrostatic and van der Waals 
contributions because P(e|</)) = P{erdw\<f>) ■ P{£eiec\4’)i where 
£ = £vdw + £eiec, i.e. the individual binding energy distribu¬ 
tions are uncorrelated. 
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